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Constant  defined  in  Equation  (26) 


Wave  speed 

Function  for  the  free  surface  given  as  z  =  F(x,y) 
Froude  number 


Ship  shape  function  given  as  y  =  f(x,z) 

Acceleration  of  gravity 
Constant  defined  in  Equation  (26) 

Wave  number 
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Water  particle  speed  on  the  free  surface  relative  to  a  ship 
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Time 
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ABSTRACT 


The  effect  on  ship  waves  of  the  nonlinear  free 
surface  boundary  condition  and  nonuniform  flow  are 
considered.  The  wave  number  is  represented  by  a 
function  of  space  and  the  wave  slope,  ak.  The  value 
of  ak  turned  out  to  be  larger  than  the  critical  value 
(0.447)  in  the  wide  area  near  the  bow,  except  when  the 
wave  crest  is  nearly  perpendicular  to  the  ship  hull; 
fortunately,  the  latter  is  the  most  significant  case. 

Yim's  previous  work  on  ray  theory  is  used  to  include 
the  full  nonlinear  effect;  the  ray  and  the  wave  phase 
of  a  Wigley  hull  are  computed  by  the  present  method 
and  are  compared  with  previous  results.  When  compared 
to  the  nonuniform  flow  effect,  the  nonlinear  free 
surface  effect  increases  further  both  the  Kelvin  wedge 
angle  near  the  bow  and  the  difference  of  wave  phase  from 
the  predictions  of  linear  theory. 

ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  was  supported  by  the  Numerical  Ship  Hydrodynamics  Pro¬ 
gram  at  the  David  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC) .  This 
program  is  jointly  sponsored  by  the  Office  of  Naval  Research  and  the  David  Taylor 
Naval  Ship  Research  and  Development  Center  under  Task  Area  RR0140302,  and  Work  Units 
1843-045  and  1542-018. 


INTRODUCTION 

Researchers  have  advanced  two-dimensional  nonlinear  wave  theory  considerably 
1* 

since  Stokes  represented  the  wave  height  by  a  series  of  products  of  wave-amplitudes 

2  3  4 

and  wave  numbers.  Lonquet-Higgins,  Schwatz,  Cokelet,  and  many  others  have 
contributed  to  developing  accurate  numerical  predictions  of  two-dimensional  wave 
forms  up  to  the  breaking  stage.  Having  assumed  that  the  wave  number  was  slowly 
varying  in  time  and  space  compared  with  the  wavelength,  Whitham'*  analyzed  three- 
dimensional  waves  using  the  two-dimensional  theory  of  nonlinear  waves.  Using 
Whitham's  idea,  Lighthill^  also  contributed  significantly  to  an  understanding  of 
three-dimensional  water  waves. 

In  a  discussion  of  a  paper  by  Gadd’7  on  ship  waves,  Lighthill  suggested  that  the 
Whitham  concept  could  be  used  to  find  a  nonlinear  correction  to  ship  waves.  When 

g 

Hogben  used  nonlinear  theory  to  determine  phase  relations  for  ship  waves,  which 
*A  complete  listing  of  references  is  given  on  page  29. 
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were  quite  different  from  the  results  of  linear  theory,  he  encountered  difficulties 
in  following  Lighthill's  advice.  Following  a  different  approach,  Hogben  applied 
two-dimensional  Fourier  expansions  with  a  large  number  of  Fourier  coefficients  to 
represent  water  waves  caused  by  a  simple  submerged  source.  The  Fourier  coefficients 
were  determined  so  as  to  satisfy  the  exact  free  surface  boundary  conditions.  Hogben 
then  obtained  wave  phases  for  various  submergences,  source  strengths,  and  Froude 
numbers. 

9 

Yim  recently  applied  a  ray  theory  of  ship  waves  and  found  the  predicted  wave 
phase  to  differ  from  the  linear  theory.  The  difference  was  caused  by  the  nonuniform 
flow  created  by  the  ship  itself.  In  ray  theory  a  ray  is  the  path  of  a  wave  energy 
packet  from  the  wave  source  such  as  the  ship  bow.  In  the  conventional  linear  ship- 
wave  theory  a  ray  is  a  straight  line  so  that,  in  relation  to  wave  propagation,  linear 
theory  neglects  the  flow  perturbation  created  by  an  advancing  ship.  However,  if  the 
flow  perturbation  is  considered  using  even  a  double  model  hull  representation,  the 
ray  is  curved,  and  some  rays  that  pass  near  the  ship  have  to  reflect  from  the  ship 
boundary.  The  Yim  analysis  considered  the  linear  free  surface  condition  as  modified 
by  the  nonuniform  flow.  The  dispersion  relation  was  linear  but  modified  to  account 
for  the  nonuniform  flow,  as  done  by  Ursell  in  1960.  A  fully  nonlinear  free-surface 
condition  includes  not  only  the  nonuniform  flow  effect  but  also  the  interaction  be¬ 
tween  wave  numbers  and  wave  amplitudes.^ 

The  present  paper  considers  both  the  effect  of  the  nonlinear  free  surface 
condition  on  ship  waves  and  the  effect  of  the  nonuniform  flow.  The  method  is  similar 
to  that  suggested  by  Lighthill.^  The  wave  number  is  represented  by  a  function  of 
space  and  the  wave  slope,  ak.  The  value  of  ak  turns  out  to  be  larger  than  the 
critical  value  (0.447)  in  a  wide  region  near  the  bow,  except  when  the  wave  crest  is 
nearly  perpendicular  to  the  ship  hull;  fortunately,  the  latter  is  the  most  signifi¬ 
cant  case.  Because  of  Yim's  previous  work  on  ray  theory,  it  was  not  difficult  to 

include  the  full  nonlinear  effect.  The  ray  and  the  wave  phase  of  the  Wigley  ship 
9 

hull  are  computed  by  the  present  method  and  are  compared  with  the  previous  results. 
For  the  first  time,  the  Whitham^  nonlinear  wave  theory  is  applied  to  the  actual  ship 
model . 
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BOUNDARY  CONDITIONS  FOR  SHIP  WAVE  PROBLEMS 
A  right-handed  rectangular  coordinate  system  is  used,  with  z  upward,  x  and  y  on 
the  mean  free  surface  with  the  origin  fixed  in  space.  At  any  time  the  fluid  is 
considered  to  be  homogeneous,  inviscid,  irrotational,  and  infinitely  deep.  When  a 
ship  advances  with  a  uniform  speed  -U  in  the  negative  x  direction,  there  exists  a 
velocity  potential  <p  that  satisfies 


V  <}> 


(1) 


with  velocity  —  V cfi ; 


-V(J)  =  ui  +  vj  +  wk  and  -V<j>  =  0  at  infinity 


The  boundary  conditions  for  <J>  on  the  free  surface,  z=F(x,y,t),  are: 


F  -  <j>  F  -  4>  F  +(<)  =  0 

t  xx  y  y  z 


(2) 


F(x,y)  <<^>yV~2V  =  0 


(3) 


w’  .e  g  is  the  acceleration  of  gravity. 

As  a  simple  approximation,  all  the  quantities  perturbed  by  the  presence  of  a  ship  in 
otherwise  uniform  flow  are  considered  sufficiently  small  so  that  the  second-order 
terms  of  Equations  (2)  and  (3)  are  negligible.  Then  Michell's  linear  boundary 
conditions  are  obtained.  When  F  is  eliminated  between  Equations  (2)  and  (3),  the 
well-known  linear  free-surface  condition  is  obtained. 


ftt  +  g4>z  =  0 


(4-1) 


If  the  coordinate  is  fixed  at  the  ship  bow,  Equation  (4-1)  becomes 


xx  ,,2  z 


4>  =  0 


(4-2) 
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The  boundary  condition  on  the  ship  surface  y  =  f(x,z)  is: 


fx 


4>  f  =  0 
z  z 


(5) 


The  Michell  thin  ship  approximation  has  been  studied  extensively  for  ship  waves  and 
wave  resistance  and  has  helped  promote  an  understanding  of  the  related  physics. 
However,  the  Michell  theory  can  predict  neither  wave  heights  nor  the  wave  resistance 
accurately,  except  when  the  ship  is  thin.  The  Michell  solution  indicates  that  the 
effects  of  beam-length  ratio  B/L  and  Froude  number  F^  =  U//gL  are  important  and 
complex,  such  that  problems  involving  high-speed  ships  should  be  treated  separately 
from  those  for  slow-speed  ships.  Also,  the  shape  of  the  ship  near  the  free  surface 
affects  the  solution  more  than  it  does  near  the  bottom  of  the  ship,  especially  for 
slow-speed  ships. 

For  a  steady  slow-ship  theory,  the  perturbation  potential  <p  and  the  free  surface 
elevation  F  are  written  as 


<t>  ”  (6) 

f  =  <;r  +  c  (7) 


2  2  2 

where  <|>r  and  =  (U  -4>rx-<l>ry) / (2g)  are  caused  by 
the  hull. 

Using  the  order  of  magnitude  assumptions: 


a  double  model 


representation  of 


<t>’  =  0(U5),  5  =  0(U4) 


and  for  cf> *  and 


/3_  3_ 
v 3x’  3y’ 


(8) 
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and  for  ^  and  ^  (f^,  |^f  |j)  =  0(1);  4>r  =  0(U),  ^  =  0(U2)  . 

The  lowest-order  free-surface  condition  for  slow  ships  to  0(U  )  has  been  derived  by 

„  .  n 

Ogilvie  as 


1  1  9  3  * 

g  ^rx(x,y’0)  9x  +<!>ry(x,y,o)  Jy  ^’(x.y.z)  +  $'2(x,y,z)  =*  D(x,y) 


on  z  =  £ 


where 


D(x,y)  =  lx  ^rx(x-yi0)^r(x’y)}  +  {<^ry(x.y>°)^r(x»y)} 


DISPERSION  RELATION 

The  Laplace  equation  is  satisfied  by  the  potential  for  an  elementary  wave  whose 
normal  to  the  crest  makes  an  angle  0  with  the  x  axis. 


kz+ia  _  kz+ik(x  cos  0+y  sin  9)-iajt 
e  =  e 


where  u  is  the  frequency  and  k  is  the  wave  number. 

Substituting  Equation  (11)  in  the  free  surface  boundary  condition.  Equation 
(4-1),  the  well-known  dispersion  relation, 

a>2  -  gk  =  0  (12) 

can  be  obtained. 

Because  the  wave  form  relative  to  a  ship  moving  at  constant  velocity  -U  is  station¬ 
ary,  the  wave  speed  c,  with  respect  to  the  calm  water  surface,  should  be  c  =  U  cos  0 
so  that 


(a  =  ck  =  lik  cos  0 


5 


This  is  a  slowly  varying  wave  number  that  satisfies  the  linear  dispersion  relation 
but  includes  the  effect  of  nonuniform  flow.  This  dispersion  relation  can  also  be 
obtained  from  the  slow-ship  free  surface  condition.  Equation  (9). 


6 


Stokes  considered  nonlinear  two-dimensional  water  waves  and  represented  the 

2 

frequency  as  a  series  in  (ak)  where  a  is  the  wave  amplitude  and  k  is  the  wave 
12 

number.  Lighthill  suggested  using  a  closed  form 

0)r2  =  gk(l+a2k2)  (18) 

as  a  good  approximation  for  a  nonlinear  dispersion  relation  because  this  is  almost 

2 

exact  for  both  critical  wave  slope  at  (ak)  =0.2  and  a  very  small  wave  slope.  The 
wave  number  may  be  obtained  from  Equations  (15)  and  (18)  and  is  given  by 

k  =  £l+a2g2/{ (U-Hij)cos  0+v^sin  0}  J g/ { (U+u^)cos  0+v^sin  0}2 

=  k^d+a2^2)  (19) 

2  2  2 
with  k^  r  g/(u  cos  +v  sin)  ,  u  =  V  +  u^,  v  =  v^  accurate  up  to  the  order  of  a  k^  . 

In  Equation  (19),  k  is  the  slowly  varying  wave  number  that  satisfies  the  non¬ 
linear  dispersion  relation  in  nonuniform  flow.  Although  the  exact  free-surface 
condition  should  be  satisfied  on  the  free  surface  and  not  on  the  mean  free  surface, 
it  is  noteworthy  that  the  wave  number  in  Equation  (19)  is  represented  by  the  double 
model  velocity  field  and  the  amplitude  relative  to  the  mean  free  surface,  z=0. 

As  can  be  seen  by  taking  the  gradient  of  the  potential  of  an  elementary  wave, 

identity  Equation  (11),  derivatives  with  respect  to  x,  y  and  z  involve  a  multipli- 

2  -2 

cation  factor  kL  =  0(gL/U  )  =  0  (F  ),  which  is  a  large  number  for  slow  ships. 

n 

Therefore,  the  derivatives  of  flow  quantities  with  respect  to  x,  y  or  z  change  the 
order  of  magnitude.  This  is  the  main  reason  for  treating  slow  ships  separately  from 
fast  ships.  Eggers  treated  the  free-surface  condition  for  fast  ships  and  derived 
a  higher-order,  free-surface  condition  on  z=0,  which  is  quite  different  from 
Equation  (9).  Then  he  obtained  a  different  dispersion  relation  which  indicated  a 
region  that  does  not  allow  waves  near  the  bow. 
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RAY  EQUATION 

When  Che  wave  number  k  is  known  as  a  function  of  space,  the  solenoidal  nature 
of  the  wave  number  vector  given  in  Equation  (14)  leads  to  the  ray  equation 

V  x  k  =  0  (20) 

where  k  =  k  cos  0  i  +  k  sin  8  j. 

Rays  for  the  wave  number  that  satisfies  the  linear  dispersion  relation,  and  includes 

9 

the  effect  of  nonuniform  flow,  have  been  investigated  in  detail  by  Yim.  The  non- 

uniform  flow  effect  curves  rays  near  the  ship.  This  is  especially  marked  when  the 

wave  crest  is  nearly  perpendicular  to  the  ship  surface.  The  ray  refracts  toward  the 

ship  and  reflects  from  the  ship  surface  to  form  the  second  caustic. 

14 

When  Adachi  analyzed  experimental  wave  height  spectra,  he  found  that  the 

measuied  amplitude  spectra  were  lower  in  value  than  predicted  by  linear  theory,  over 

the  whole  range  of  0,  and  much  lower  at  values  of  (i  less  than  20  degrees.  The 
sheltering  effect^"*  has  been  shown  to  reduce  wave  resistance  15-25%  below  that  pre¬ 
dicted  by  linear  theory,  which  is  consistent  with  reduced  wave  amplitudes.  However, 
it  was  not  known  why  the  amplitudes  were  particularly  reduced  at  values  of  0  less 

than  20  degrees.  If  the  wave  reflection  of  elementary  waves  for  20  degrees  is 

considered,  the  experimental  phenomena  can  be  understood.  Although  the  flow  field 
near  the  ship  may  be  further  illuminated  experimentally  by  comparing  a  long  wedge 

type  of  ship  bow,  which  does  not  allow  ray  reflection,  and  a  rounded  bow  which  in- 

.  9 

duces  multiple  ray  reflections.  A  bow  bulb  reduces  the  nonuniform  flow  and  also 
prevents  reflection  if  a  proper  size  bulb  is  used.^ 

When  only  the  conventional  linear  dispersion  relation  in  regard  to  wave  number 
is  considered  without  the  nonuniform  flow  effect,  all  the  rays  will  be  straight  1 ines , 
and  the  ray  theory  adds  no  useful  information  to  that  obtained  bv  applying  the 
linear  wave  theory.  To  obtain  useful  information,  researchers  must  consider  non- 
uniform  flow  and/or  a  nonlinear  dispersion  relation.  Although  ray  theory  with  a 
nonlinear  dispersion  relation  has  never  been  considered  for  ray  paths,  it  may  not  be 
too  difficult  to  include  the  effects  of  both  nonlinear  dispersion  and  nonunitorm 
f  low. 
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NONLINEAR  RAYS 

From  Equation  (20) 

3kL  3k2 

3y  3x 

where  =  k  cos  0,  k2  =  k  sin  0 
or 


(21) 


k  cos  0  -k  sin  0  -k  sin  00  -k  cos  00  =0 

y  x  y  x 


(22) 


Differentiation  of  Equation  (19)  yields 


k  =  (l+3a2k2)  k.  +  2  a  k  3  a 

x  8.x  SI  x 


k  =  (l+3a2k2)  k.  +  2  a  k  3  a 
y  Jly  l  y 


(23) 


and  substituting  k^  and  k^  from  Equation  (19)  into  Equation  (22),  we  have 
{ (l+a2k^2)sin  8(u  cos  0+v  sin  0)+2(l+3a2k^2)cos  0(-u  sin  0+v  sin  0) 

O  O  O 

-2ak^  a^  cos  0(u  cos  9+v  sin  0)}  +  { (1+a  k^  )  cos  0(u  cos  0+v  sin  0) 

2  2  2  3  C) 

-2(l+3a  )  sin  0(-u  sin  0+v  cos  0)+2ak^  sin  0(u  cos  0+v  sin  0);  ^ 


=  2(l+3a2k^2)  sin  0  |cos  0  -|^  +  sin  9 


-2(l+3a2k^2)  cos  0  |cos  0  +  sin  0  -|^j  -2ak^2  (sin  0  -  cos  0  j 


3y  / 


•  (u  cos  0+v  sin  0) 
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This  first  order  partial  differential  equation  for  0  is  equivalent  to  the  simulta¬ 
neous  ordinary  differential  equations 


~  ~  (l+-a2k^2)  cos  0(u  cos  0+v  sin  0)  -2  ( l+-3a2k^2)  sin  6(-u  sin  0+v  cos  9) 


+-  2  a  k„  a„  sin  0(u  cos  0+v  sin  0) 

Ji,  O 


&  = 
dt 


2  2 

( 1+a  k^) 


sin  0(u  cos  0+v  sin  6) 


2  2  2 
+  2(l+-3a  k^  )  cos  0(-u  sin  0+v  cos  0)  -2  a  k^  cos  0(u  cos  0+-v  sin  9) 

(24) 

d6  t  f,  .  o  2.  2.  .  0  /  A  9u  .  A  9v  \ 

dF  =  2(1+3a  k£  )  sin  0  lcos  9  ^  +  Sln  0  9^) 

-  2(l+-3a2k^2)  cos  0  ^cos  0  +  sin  8  j 

-  2  a  k^2  ^sin  0  ^  -  cos  0  j  (u  cos  8+v  sin  0) 


with  initial  conditions  placed  on  x,  y,  and  0. 

At  present,  there  is  no  other  choice  of  wave  amplitude,  a,  except  that  of  linear 

theory  or  that  of  the  linear  theory  modified  by  the  sheltering  effect.23  Because  in 

2 

the  ray  equation  the  amplitude,  a,  appears  only  as  a  ,  the  linear  amplitude  may  be 

acceptable,  except  near  the  first  caustic  at  the  Kelvin  wedge,  where  a  is  singular. 

The  linear  wave  number  kp  is  very  large  near  the  ship  bow  where  the  stagnation  point 

t  2 

is  located.  Thus,  the  region  where  (ak)  >0.2  will  be  near  the  bow  and  near  the 

2 

first  caustic.  Because  there  is  no  two-dimensional  wave  with  values  of  (ak)  larger 

than  the  critical  value  of  0.2  with  the  present  theory,  it  is  impossible  to  consider 

2 

nonlinear  ship  waves  with  (ak)  >  0.2. 

The  linear  wave  amplitude  far  from  the  wave  source  can  be  obtained  by  the  method 
of  stationary  phase,22  where  the  amplitude  is  a  function  of  0  and  proportional  to 
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r  ;  r  represents  the  distance  between  the  field  point  and  the  wave  source.  Because 
the  wave  action  contained  in  a  ray  tube  is  unchanged,  when  the  far  field  wave  energy 
and  the  ray  paths  are  known  the  wave  energy  can  be  obtained.  Wave  action  is  the 
wave  energy  divided  by  the  relative  frequency,  and  the  wave  energy  is  a  function  of 
wave  amplitude. 

The  wave  phase  s  can  be  obtained  in  the  same  way  as  for  only  the  nonuniform  flow 
effect  by 


ds  =  k  •  dr  =  k„ 


(1+a^k^)  cos  8  dx  +  k^ 


(l+a^k.^)  sin  8  dy 


(25) 


However,  when  (ak)  >  0.2,  the  wave  breaks  and  the  physics  of  the  wave  propagation 

and  the  energy  relation  may  change.  Nevertheless,  the  maximum  phase  change  can  be 

2 

obtained  by  simply  using  a  value  of  0.2  for  (ak)  in  Equation  (25),  which  yields  1.2 
times  the  value  of  the  phase  obtained  when  only  the  nonuniform  flow  is  considered. 

When  the  linear  wave  amplitude  is  used  to  solve  Equation  (24),  the  choice  of 
initial  values  of  (x,y,0)  poses  a  serious  problem.  Previous  experience  has  shown 
that,  given  the  nonuniform  flow  effect,  many  rays  are  outside  of  the  first  caustic 
predicted  by  linear  theory,  where  no  linear  wave  is  supposed  to  exist.  In  addition, 
the  value  of  linear  wave  amplitude  is  infinite  at  the  first  caustic.  In  both  linear 
and  nonlinear  theories  the  wave  amplitudes  at  x  for  the  same  value  of  9^  may  be 
assumed  to  be  about  the  same  for  the  present  approximation.  However,  in  the  present 
nonlinear  theory  rays  at  x  -*■  00  can  be  obtained  only  when  the  approximate  amplitude 
is  known  all  the  way  to  x  +  °"  from  the  initial  point  near  the  origin.  Therefore, 
rays  due  to  nonuniform  flow  only,  without  amplitude  interaction,  were  computed  first 
together  with  the  value  of  6  . 

The  0^  value  was  then  used  to  compute  linear  wave  amplitudes  by  the  stationary 
phase  method.  Because  rays  are  assumed  to  start  from  a  point  wave  source,  their 
amplitudes,  as  obtained  by  either  the  energy  method  or  the  stationary  phase  method, 
are  singular  at  the  initial  point  or  at  the  wave  source.  The  stationary  phase 
method  is  used  for  the  asymptotic  value  of  the  amplitude  at  large  x  values.  How¬ 
ever,  if  rays  are  considered  to  be  straight  lines,  the  same  value  of  amplitude  near 
the  wave  source  results  from  either  the  stationary  phase  method  or  the  energy  method 
for  a  constant  wave  number.  Therefore,  in  the  present  analysis  the  stationary  phase 
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method  was  applied  in  all  cases  to  obtain  the  amplitudes  and  their  derivatives  for 
use  in  computing  nonlinear  ray  paths  and  0.  The  final  0  value  at  x  -*•  00  for  the  non¬ 
linear  ray  paths  may  be  different  from  0^  which  is  used  for  the  linear  wave  ampli¬ 
tude.  However,  the  linear  amplitude  changes  slowly  at  points  not  too  near  the  wave 
source  for  small  0^  values  where  the  present  calculation  of  nonlinear  ray  paths  is 
meaningful.  Therefore,  errors  of  higher  order  may  be  induced  by  using  the  approxi¬ 
mate  values  of  wave  amplitude.  Some  numerical  results  given  below  for  a  Wigley  hull 
seem  to  indicate  that  this  conjecture  is  valid. 

EXAMPLES  AND  DISCUSSION 

Using  Equation  (24),  nonlinear  rays  are  computed  for  a  Wigley  hull  whose  source 
distribution  is  represented  by 


m  =  b^  (l-2x)  j-1-  j  on  0  <  x  <  1,  -  h<z<h  (26) 

Appendix  A  contains  equations  for  the  wave  amplitudes  and  their  derivatives  for  the 

Wigley  hull.  Equations  for  the  flow  field  to  be  used  in  the  ray  equation  are  given 

in  Appendix  B.  These  equations  are  exactly  the  same  as  those  in  Reference  9  when 

several  misprints  in  Reference  9  have  been  corrected. 

In  Figures  1  through  3,  ray  paths  computed  according  to  Equation  (24)  are 

9 

compared  with  earlier  results  derived  by  Yim  in  which  only  the  nonuniform  flow 
effect  was  considered  without  the  effect  of  amplitude  interaction.  Inclusion  of  the 
amplitude  interaction  seems  to  push  the  ray  away  from  the  ship  when  0  | 0  [  <  35  deg 
(0.61  radians)  just  as  the  nonuniform  flow  effect  did  to  the  conventional  linear  ray 
theory.  When  |0  |  >  35  deg  (0.61  radians),  the  amplitude  interaction  effect  pulls 
the  ray  toward  the  ship. 

In  Figure  3  the  ray  paths  for  0^  =  -0.307  and  -0.292  radians  are  almost 

identical,  although  the  linear  amplitudes  are  used  for  0^  =  -0.246  and  -0.0909 

radians,  respectively.  This  means  that  the  linear  amplitudes  for  computation  of  the 

nonlinear  ray  theory  are  reasonably  insensitive  to  0^  for  small  6^  values. 

9 

In  Figures  4  through  6,  82(6)  denotes  the  radian  phase  difference  divided  by 
k^  of  waves  as  predicted  by  the  ray  theory  and  linear  theory.  At  far  down  stream 
the  wave  number  is  constant,  and  for  0-0  the  wave  number  equals  k.  Therefore, 
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Figure  4  -  Phase  Difference,  s^(0),  between  Linear  and  Nonlinear 
Theory  of  the  Wigley  Hull  (h=0.2,  h=0.0625,  F^=0.25)  at  x=2 


Figure  5  -  Phase  Difference,  s^fO),  between  Linear  and  Nonlinear 
Theory  of  the  Wigley  Hull  (b=0.2,  h=0.03,  Fn=0.25)  at  x=2 
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0(x*2t 

Figure  6  -  Phase  Difference,  82(6),  between  Linear  and  Nonlinear 

Theory  of  the  Wigley  Hull  (b=0. 2,  h=0.0625,  F  =0.28  and  0.25) 

at  x=0. 5 

s2(0)  may  be  considered  to  be  the  fractional  wavelength  nondimensional ized  by  the 

ship  length  converted  from  the  wave  phase  difference  because  the  wavelength  is 

2ir/k  .  As  Figures  4-6  show,  the  effect  of  amplitude  interaction  on  the  wave  phase 
o 

is  considerably  larger  than  the  effect  of  nonuniform  flow.  However,  the  total 

magnitude  of  the  phase  differences  from  the  linear  theory  seems  to  be  quite  close  to 

1  ft 

experimental  results1  showing  phase  differences  of  3-5%  of  the  ship's  length  (L) 
along  the  ship  hull. 

Q 

Hogben  calculated  the  phase  difference  between  the  nonlinear  and  linear  waves 
caused  by  a  submerged  point  source.  Because  the  case  involved  a  submerged  body,  he 
was  able  to  demonstrate  the  phase  advancement  in  the  nonlinear  theory  without  having 
to  consider  the  sheltering  effect.  In  the  present  theory  the  phase  advancement  was 
clarified  further  by  separating  it  into  contributions  caused  by  amplitude  interaction 
and  by  nonuniform  flow.  In  Figure  6,  the  values  of  phase  differences  normalized  by 
the  wave  number  change  little  in  the  Froude  number  range  0.20  <  Fn  <  0.27  although, 
as  Hogben  shows,  these  phase  differences  are  slightly  less  for  higher  Froude  numbers. 

In  the  nonlinear  theory  the  principle  of  superposition  does  not  work,  as  is  well 
known.  The  present  results,  however,  derive  from  a  nonlinear  analysis  of  each 
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elementary  wave.  The  elementary  waves  form  the  bow  or  stern  wave  by  superposition. 
When  the  regular  wave  integral  is  evaluated  bv  the  stationary  phase  method,  the  wave 
usually  is  taken  to  be  the  superposition  of  two  wave  systems,  transverse  and 
divergent.  However,  because  the  divergent  waves  are  known  to  be  of  almost  zero 
amplitude  near  the  ship  hull,  only  transverse  waves  are  significant.  In  this  sense, 
the  present  evaluation  of  rays  near  the  hull  surface  is  meaningful.  At  some  distance 
from  the  hull,  the  predicted  approximate  behavior  of  ray  paths  and  phase  differences 
help  one  to  understand  discrepancies  between  the  predictions  of  linear  theory  and 
the  experimental  results. 

The  effect  of  nonlinear  amplitude  interaction  on  reflecting  rays  has  been  tested 
numerically  as  shown  in  Figure  7.  The  amplitude  interaction  cannot  obstruct  the 
reflection  of  rays  from  the  ship  surface.  Reflecting  rays  are  pushed  outboard  due 
to  the  amplitude  interaction,  causing  the  shape  of  the  second  caustic  to  change 
considerably .  The  paths  of  rays  and  reflected  rays  are  highly  dependent  on  the  hull 
form  through  the  nonuniform  flow  and  wave  amplitude  interaction  effects.  In 
addition,  near  the  caustic  even  the  wave- induced  flow  field  interacts  with  the  wave 
propagation.  Therefore,  it  may  be  quite  difficult  to  verify  quantitatively  the 
location  of  the  second  caustic  observed  in  the  towing  tank.  Nevertheless,  the  in¬ 
fluence  of  the  wave-induced  velocity  near  the  second  caustic  can  be  considered  in  the 
following  way. 


0  0.1  0.2  0.3  0.4  0.5 


X 

Figure  7  -  Amplitude  Interaction  on  Ray  Reflections  from  the  Wigley  Hull 

(b=0. 2,  h=0. 0625,  Fn=0.2) 
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The  ray  equations  which  include  the  amplitude  interaction  can  be  obtained  from 
Equation  (24)  by  neglecting  the  quadratic  amplitude  terms: 


d6  . 

dt  =  2 


(  a  •  a  o(  a  ,  .  a  9v  \ 

^cos  0  4*  sin  0  J  -  2^cos  0  +  sin  0 


d  x 

—  =  2  sin  0  (u  sin  9-v  cos  0)  +  cos  0  (u  cos  6+v  sin  6) 
d  t 


The  first  equation  can  be  written  in  a  vector  form 


d6 

dt 


(27) 


where  9/9t^  is  along  the  tangent  to  the  crest  of  the  wave,  and  the  vectors  are  given 
by 


t^  =  i  sin  9  -  j  cos  0 
n  =  i  cos  0  +  j  sin  8 


and 


q  =  ui  +  vj 

9 

Because  the  reflecting  wave  crest  is  almost  perpendicular  to  the  ship  hull  and  q 
near  the  ship  is  almost  parallel  to  the  ship  hull,  n  and  q  are  almost  parallel,  i.e., 
n  •  q  '  /q/.  If  the  wave  height  increases,  / q /  will  decrease  according  to  linear 
wave  theory.  Therefore,  if  the  ray  crosses  the  second  caustic  into  the  large 
amplitude  wave  region,  d0/dt  will  be  negative  and  0  will  decrease  according  to 
Equation  (27).  Conversely,  if  the  reflected  ray  crosses  the  second  caustic  out  of 
the  large  amplitude  wave  region,  dO/dt  will  be  positive  and  0  will  increase.  Then 
the  reflected  ray  will  go  farther  from  the  ship  hull.  This  will  cause  the  second 
caustic  to  occur  farther  from  the  ship  hull,  thus  leaving  a  wider  region  of  breaking 
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waves  near  the  ship.  Therefore,  the  actual  breaking  wave  region  will  be  larger  than 

that  predicted  by  the  nonuniform  flow  effect  only.  The  energy  flux  of  the  elementary 

bow  waves  near  0^  =  0  that  can  reflect  from  a  ship  hull  is  limited  to  a  certain  range 

of  the  initial  value  of  0.  The  range  depends  upon  the  bow  shape.  If  the  bow  shape 

9 

is  a  wedge,  there  will  be  no  reflection.  If  it  is  a  large  round  bow,  there  may  be 
a  large  range  of  0  whose  rays  will  reflect  from  the  ship  hull.  However,  near  the 
stagnation  point  and  the  boundary  layer  the  flow  does  not  vary  slowly  enough  for  the 
present  theory  to  be  applied.  Therefore,  ships  of  large  beam  length  ratio  have  a 
large  curvature  near  the  bow  and  show  in  the  towing  tank  a  few  of  the  distinctive 
phenomena  predicted  here.  However,  for  a  thin  model  in  a  towing  tank,  the  reflected 
wave  energy  may  be  too  weak,  and  the  second  caustic  be  too  close  to  the  ship  hull  to 
overcome  the  boundary  layer  effect.  In  this  case  linear  theory  may  predict  rather 
well.  To  conclude,  the  nonlinear  effects  of  amplitude  interaction  on  ray  paths  and 
phases  of  ship  waves  have  been  investigated  using  the  Stokes*  nonlinear  wave  formula 
and  the  Whitham'’  nonlinear  wave  theory  in  conjunction  with  ray  theory.  For  the  ray 
path  the  nonuniform  flow  effect  is  greater  than  the  effect  of  the  amplitude  inter¬ 
action,  but  for  the  phase  the  latter  effect  is  greater  than  the  former. 
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APPENDIX  A 

LINEAR  WAVE  AMPLITUDE 


In  the  ray  equation  the  computed  value  of  the  linear  wave  amplitude  a  is  used. 
The  height  of  the  regular  bow  wave  caused  by  the  source  distribution  m  for  the 
Wigley  hull  is 


C,  =  4 


tt/2  0  1 

J  1/ 

-tt/2  -h  0 


k  z . 
o  1 


m(x1,z1)  e 


2  2 

sec  0  +  i  sec  0{(x-x^)eos  0+y  sin  0) 

k  sec"*  0  dx.  dz  d0 

o  11 


where 


then 


k  =  gL/U  =  F 
oa  n 


m(x^,z^)  =  b-^l-x.^  ) 


tt/2 


5  =  4  b 


1  f  E1  (  S"  +1)  e 

J  \k  sec  8  / 

-tt/2  '  o  7 


i  k  (x  sec  0+y  sec  0  sin  0) 
o 

d0 


E  =  -f - 


2  +^ 


-k  h  sec  0 
o 


fl+ 


'1  k  v2  ,  3  4.  ,.2  2,1..  2fi 

o  h  k  sec  0  h  k  sec  0  \  k  h  sec  0 

O  O  NO' 


Using  the  method  of  stationary  phase,  the  amplitude  can  be  obtained  as 


a  =  4  b2  Ej 


,2  2  , 
ik  sec  8 
N  o 


i+1)(kor|F”»l)l) 


1/2 
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where 


x  =  r  cos  6,  y  =  r  sin  & 


F(6)  =  cos  6  sec  8  +  sin  6  sec  9  sin  8 


0^  satisfies  the  equation 


2 


2 

tan 


+  cot  6  tan  0^+1 


0 


or 


tan  6  = 


tan  8^/(2  tan  0^+1) 


and 


F" (0L)  =  sec  0l(1-2  tan201)  {(2  tan  9^1)+  tan  0^ 


-1/2 


where 


F" ' (0) 


1/2 
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8  tan  9 

2  3  A  a 

h  k  sec  o 
o 


-k  h  sec26 
o 

b  e  tan  8 
k 

o 


-k  h  sec  8 

A  e  _ _ _  cos^  0  tan  0 

k3h2 

o 


ax  a9  3x  9x 


00  ,  3a 


00  0a 

ay  =  a0  0y  + 


where  ooly  r  U  coostdered  to  be  .  faction  of  x  or  y  t«  ^  »-  W3y' 
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APPENDIX  B 

FLOW  FIELD  OF  A  WIGLEY  HULL 


For  the  computation  of  ray  paths  of  a  Wigley  hull,  the  flow  velocity  and  its 
derivatives  u,  v,  u^,  v  ,  u  are  needed  on  (x,y,o).  A  Wigley  hull  has  the  double 
model  source  distribution 

m  =  b^ (-2x^+1) 
in 


0  <  x!  <  1*  yl  =  °»  h  >  zi  >  -  h 


Thus , 


-h  -1 

-  u  (x,y,o)  -  2  1  J  m  (i)  dx1  dZl  -  1 
0  0 


h  1 


=  -  2  f  J  m  dxl  dzl  - 

Jo  *0  1 


=  2  b, 


.h  1-  -=x  -hi 


J  7(7^Tydzl+J  rCx^oj-  dzl 


-4  b, 


2 

3  h  log  {x^-x+r(z 


rh  (x  -x)  y2 

i=h>}-  7  V  2:  ~ 

0  (zi  +y  }  r 


dz. 


fh  (Xj-x)  1  fh  4  (x2~x)  (Xj-xXz^-y2) 

'J  r  dzl  '  772  J  y  .  2  2,  r  dzl 

0  3  h  0  (Z1  +y  >  r  J 


25 


aacsmm  ho*  bumm<q *  num 
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OJ  1 1— • 


h  1 


3 2  a 


Uy  =  2  /  J  “  9y3x 

0  0 


(r )  dxl  dzl 


-  2  bxy 


+  +  1 


{(l-x)2+y2)r(l,h)  r(1-h)h  r(°'h)h 
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h  (1-x) 


Q  (z12+y2)r(x1=l) 
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9  ?  d^  1 
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h  r 

1  f  1-x  ,  1  (  x  , 

2  J  r(x  =1)  dZl  '  ,2  J  r(x  =0)  d  1 

“  rv  *■  “  n 


-bi  (6 


Vy  -  ^6y  ■  h2 


2+  ^ 


i  i  !  (xrx)2-y2 

i(1-2x)U^  +  ~; — rr 

1  1  (x^-x)y 


-x)2+y2}r(z^=h) 


o  (zi2+y  )r  3  {(X1 

(l-2x) (  - j — 2 - 

\  (x1-x){(x-x1)  +y  }r(zx 


hrU^h) 


(h2+y2) (x^-x)y2/ 


-6  b^y 


xr° 


i/Jl 

3h2  \{  (Xj^ 


2yz+(l-2x) (x1~x)h 


3h“  \( (x,-x)2+y2)r(z1=h) 
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+  2  log 


+ 


y 


h+r (z^=h) 
r(z1=0) 


2h 


r (z1=h) 


2(x1~x) (l-2x) 
3h2 


f  d2l 

{  (z. 2+y2) r 


1  1 


v(x. 


x^=0 


where 


r(a,b)  =  r(x1=a,  z^b) 


In  these  expressions  the  integrals 


are  known  in  closed  form. 
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